Quantum vortex reconnections

We study reconnections of quantum vortices by numerically solving the governing Gross-Pitaevskii equation. We find that the minimum distance between vortices scales differently with time before and after the vortex reconnection. We also compute vortex reconnections using the Biot-Savart law for vortex filaments of infinitesimal thickness, and find that, in this model, reconnection are time-symmetric. We argue that the likely cause of the difference between the Gross-Pitaevskii model and the Biot-Savart model is the intense rarefaction wave which is radiated away from a Gross-Pitaeveskii reconnection. Finally we compare our results to experimental observations in superfluid helium, and discuss the different length scales probed by the two models and by experiments.Comment: 23 Pages, 12 Figure

Interpolating discrete advection-diffusion propagators at Leja sequences

We propose and analyze the ReLPM (Real Leja Points Method) for evaluating the propagator phi(DeltatB)nu via matrix interpolation polynomials at spectral Leja sequences. Here B is the large, sparse, nonsymmetric matrix arising from stable 2D or 3D finite-difference discretization of linear advection-diffusion equations, and phi(z) is the entire function phi(z) = (e(z) - 1)/z. The corresponding stiff differential system y(t) = By(t) + g,y(0) =y(0), is solved by the exact time marching scheme y(i+1) = y(i) + Deltat(i)phi(Deltat(i)B)(By(i) + g), i = 0, 1,..., where the time-step is controlled simply via the variation percentage of the solution, and can be large. Numerical tests show substantial speed-ups (up to one order of magnitude) with respect to a classical variable step-size Crank-Nicolson solve

Vortex reconnections in atomic condensates at finite temperature

The study of vortex reconnections is an essential ingredient of understanding superfluid turbulence, a phenomenon recently also reported in trapped atomic Bose-Einstein condensates. In this work we show that, despite the established dependence of vortex motion on temperature in such systems, vortex reconnections are actually temperature independent on the typical length/time scales of atomic condensates. Our work is based on a dissipative Gross-Pitaevskii equation for the condensate, coupled to a semiclassical Boltzmann equation for the thermal cloud (the Zaremba-Nikuni-Griffin formalism). Comparison to vortex reconnections in homogeneous condensates further show reconnections to be insensitive to the inhomogeneity in the background density.Comment: 6 pages, 4 figure

High-order time-splitting Hermite and Fourier spectral methods

In this paper, we are concerned with the numerical solution of the time-dependent Gross-Pitaevskii Equation (GPE) involving a quasi-harmonic potential. Primarily, we consider discretisations that are based on spectral methods in space and higher-order exponential operator splitting methods in time. The resulting methods are favourable in view of accuracy and efficiency; moreover, geometric properties of the equation such as particle number and energy conservation are well captured. Regarding the spatial discretisation of the GPE, we consider two approaches. In the unbounded domain, we employ a spectral decomposition of the solution into Hermite basis functions: on the other hand. restricting the equation to a sufficiently large bounded domain, Fourier techniques are applicable. For the time integration of the GPE, we study various exponential operator splitting methods of convergence orders two, four, and six. Our main objective is to provide accuracy and efficiency comparisons of exponential operator splitting Fourier and Hermite pseudospectral methods for the time evolution of the GPE. Furthermore, we illustrate the effectiveness of higher-order time-splitting methods compared to standard integrators in a long-term integration

Bivariate polynomial interpolation on the square at new nodal sets

As known, the problem of choosing ``good'' nodes is a central one in polynomial interpolation. While the problem is essentially solved in one dimension (all good nodal sequences are asymptotically equidistributed with respect to the arc-cosine metric), in several variables it still represents a substantially open question. In this work we consider new nodal sets for bivariate polynomial interpolation on the square. First, we consider fast Leja points for tensor-product interpolation. On the other hand, for classical polynomial interpolation on the square we experiment four families of points which are (asymptotically) equidistributed with respect to the Dubiner metric, which extends to higher dimension the arc-cosine metric. One of them, nicknamed Padua points, gives numerically a Lebesgue constant growing like log square of the degree

A minimisation approach for computing the ground state of Gross\u2013Pitaevskii systems

In this paper, we present a minimisation method for computing the ground stateof systems of coupled Gross\u2013Pitaevskii equations. Our approach relies on a spectral decomposition of the solution into Hermite basis functions. Inserting the spectral representation into the energy functional yields a constrained nonlinear minimisation problem for the coefficients. For its numerical solution, we employ a Newton-like method with an approximate line-search strategy. We analyse this method and prove global convergence. Appropriate starting values for the minimisation process are determined by a standard continuation strategy. Numerical examples with two and three-component two-dimensional condensates are included. These experiments demonstrate the reliability of our method and nicely illustrate the effect of phase segregation

Qualidade de cocção de grãos de arroz translúcidos e gessados.

objetivou-se avaliar o comportamento de cocção dos grãos de arroz translúcidos e gessados e seus aspectos de qualidade

Características físicas dos grãos de arroz translúcidos e gessados.

Objetivou-se neste trabalho avaliar as características de transparência, brancura e grau de polimento dos grãos de arroz translúcidos e gessados

A massively parallel exponential integrator for advection-diffusion models

This work considers the Real Leja Points Method (ReLPM) for the exponential integration of large-scale sparse systems of ODEs, generated by Finite Element or Finite Difference discretizations of 3-D advection-diffusion models. We present an efficient parallel implementation of ReLPM for polynomial interpolation of the matrix exponential propagators. A scalability analysis of the most important computational kernel inside the code, the parallel sparse matrix\u2013vector product, has been performed, as well as an experimental study of the communication overhead. As a result of this study an optimized parallel sparse matrix\u2013vector product routine has been implemented. The resulting code shows good scaling behavior even when using more than one thousand processors. The numerical results presented on a number of very large test cases gives experimental evidence that ReLPM is a reliable and efficient tool for the simulation of complex hydrodynamic processes on parallel architectures

Dichotomy between urban and rural areas: statistical data may not reveal the synergy between these two existing spaces.

An analysis of several indicators, socio-economic and environmental, through the Dashboard of Sustainability is possible to tell which category has the highest rate of farmer sustainability, whether smallholders or monoculture. However, the secondary data available in Brazil today does not support a thorough analysis of the participation of each actor and to which the interconnection between the actors and their synergy in local economic activity. Since, given the narrowing between urban and rural, the statistical data available are not able to demonstrate the extent to which gives the rural-urban dichotomy. Thus, it is relevant to point out and discuss ways to provide consistent statistical data and be, in fact, able to demonstrate the local reality of a region within the welfare actors.ICAS 2013